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Abstract 

Energy levels of the double-A hypernuclei aa^s, aaLi, jSj-'h aa^i, aa-^^ Ajfi^ predicted 
on the basis of the a+x+A+A four-body model with x = n,p, d, t,^He and a, respectively. Interac- 
tions between the constituent particles are determined so as to reproduce reasonably the observed 
low-energy properties of the a + x nuclei (^He, ^Li, ^Li, ^Li, ^Be, ®Be) and the existing data of A- 
binding energies of the X -I- A and a -I- a; -I- A systems (aH, ^H, %lle, aHc, aLI, a^I, i^i' A-^^' A-^^)- 
Here, an effective AA interaction is constructed so as to to reproduce, within the a -|- A -|- A model, 
the Baa of aa-^^ which was discovered recently in the NAGARA event of the emulsion experi- 
ment. With no adjustable parameters for the a + ,t + A + A systems, .Baa of the ground and 
bound excited states of the double-A hypernuclei with A = 7 — 10 are accurately calculated with 
the Gaussian-basis coupled-rearrangement-channel method. The Demachi- Yanagi event, observed 
recently for aa^6, is interpreted as observation of its 2+ excited state on the basis of the present 
calculation. Structure change of the a + x core nuclei due to the participation of the A particles is 
found to be substantially large and it plays an important role in estimating the AA bond energies 
of those hypernuclei. 
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I. INTRODUCTION 



A recent finding of the double-A liypernucleus aaHc, wliicli is called as NAGARA event 
in the KEK-E373 experiment [|I|, has a great impact not only on the study of baryon-baryon 
interactions in the strangeness S = —2 sector but also on the study of dynamics of many- 
body systems with multi-strangeness. The importance of this event is attributed to the 
well-defined explanation of the process and the high quality experimental value of the AA 
binding energy Saa = 7.25 ± 0.19±oJf MeV [0, which leads to a smaller AA binding, 
Ai?AA = 1-01 ± 0.20±Q}f MeV, than the previous understanding. Sometimes the emulsion 
events include ambiguities related with serious difficulty of identifying emission of neutral 
particles such as neutrons and 7-rays. In the NAGARA event, however, the production of 
aaHc has been uniquely identified free from such an ambiguity on the basis of the observation 
of sequential weak decays. 

Historically, in the 1960 's, there appeared two reports on the observation of double- 
A hypernuclei, AA^e and aaHc 0, but the reality of the latter case was considered 
doubtful Two decades later the modern emulsion-counter hybrid technique has been 
applied in the KEK-E176 experiment [Q, in which a new double-A hypernucleus event was 
found but no unique identification was given so far: One explanation as l\Be leads to a 
repulsive AA interaction (ASaa < 0), while the other possibility involving l\B leads to an 
attractive AA interaction @, 0. If the latter is the case, the extracted strength of the AA 
interaction is attractive with AB\\ ~ 4 MeV. Although the latter option seems consistent 
with the old data of l^Be [^, the substantially attractive AA interaction has not been 
convincing. 

In the strangeness nuclear physics, the most fundamental problem is to recognize various 
facets of interactions among octet baryons (A^, A, S, H) in a unified way. Our detailed 
knowledge for the S = NN sector is based on the rich data of A^A^ scatterings as well as 
nuclear phenomena. Recent studies for S = —1 many-body systems such as A hypernuclei 
have clarified interesting features of AA^ and J^N interactions in spite of scarce data of 
the free-space scatterings. On the other hand, for baryon-baryon interactions with S = 
—2 sectors, concerned presently, experimental information has been highly limited due to 
the extreme difficulties of two-body scattering experiments. Therefore the observed AA 
bond energies of double-A hypernuclei should be the most reliable source for the S = —2 
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interaction, and such data play a decisive role in determining the strength of underlying AA 
interactions. 

In view of this relevance and the experimental situation, the NAGARA event is an epoch- 
making one which provides us with a new and firm basis for understanding the double-A 
hypernuclei. In recent years several experiments to produce S = —2 systems (El 76 and 
E373 at KEK, E885 and E906 at BNL) have been performed and some of the data analyses 
are still in progress to get novel information on the S = —2 interactions. 

In this exciting situation of the experimental study, it is needed to perform careful the- 
oretical calculations of double-A hypernuclei with refreshed viewpoints. As one of the mo- 
tivations of the present work, we think it necessary and timely to put the NAGARA data 
of aaHg binding energy as a new standard basis for a systematic study of a series of several 
double-A species. Secondly, in order to extract information on the AA interactions precisely, 
here we emphasize that hypernuclear calculations should be complete and realistic enough 
to leave structural ambiguity as negligibly as possible. All the dynamical changes due to 
successive A participation should be also taken faithfully. To meet these requirement we ex- 
plore light p-shell double-A hypernuclei (A = 6 — 10) comprehensively using the microscopic 
three- and four-body models. Thirdly, by these systematic and realistic calculations, we 
will give reliable prediction of not only the ground-state binding energies but also possible 
excited-state energies, which encourages double-A hypernuclear spectroscopic study in near 
future. 

So far several cluster models have appeared to estimate the ground-state binding energies 
of double-A species: Based on the old data of AA^e and aa^g, Takaki et al. P| applied a 
simplified version of the a + x + A + A cluster model to A = 6 — 10 systems in which they 
put several angular momentum restriction and neglected rearrangement channels. Bodmer 
et al. IP, performed variational Monte Carlo calculations for a + A + A and a + a + A + A 
to investigate consistency between the AA-binding energies, -BaaIaaHs) and i^AAlAA^e), 
although their old data should be now updated. In the latest stage of this work, we encoun- 



tered with the Faddeev-Yakubovsky calculations of aa-^g and aa-^g by Filikhin and Gal |Tl 



who restricted the equations within the s-wave. They compared the results with our pre- 



vious cluster-model calculation |T3[ which was performed with wider model space but the 
stronger AA interaction strength. In our previous work [|T3] , AA binding energies have been 
calculated for aaHg and aa-^g in the framework of the a -|- A + A three-body model and the 
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a + a + A + A four-body model, respectively, where the adopted AA interaction is taken 
to be considerably attractive on the basis of the traditional interpretation for the double-A 
events. 

In the present work, by noting the importance of the NAGARA data, we extend this four- 
body model to more general cases consisting ofa + x + A + A systems with x = n,p, d, t,^He 
and a (aaHg, aa^I, aaLI, aa^I, aa^g and AA^e), where nuclear core parts are quite well 
represented by a + a; cluster models (for example in Ref. fl^)- Here we remark that the 
extensive calculations are presented for the first time for A = 7 — 9 double-A species and 
that the old predictions for aaHc and aa^g have been updated in a unified way. The four- 
body calculations are accurately performed by using the Jacobian-coordinate Gaussian-basis 



method of Refs. |]15|, |T6|, |T^, |T§, 0, ^ with all the rearrangement channels taken into 
account. In our model, structure changes of nuclear cores caused by added one and two A 
particles are treated precisely. Namely, we take into account the rearrangement effects on 
AA bond energies induced by changes of nuclear cores. It is worthwhile to point out that the 
important effects of core-excitations and core-rearrangement are lacking in the frozen-core 
approximation used often for calculations of double-A hypernuclei. 

In our model, it is possible to determine the ax and Ax interactions so as to reproduce 
all the existing binding energies of subsystems (a + x, x + A, a + x + A and a + A + A) in an 
a + x + A + A system, where that of a + A + A is given by the NAGARA event. This feature 
is important to discuss the energy levels of the double-A hypernuclei and to extract the AA 
interactions because the ambiguities of A^A^ and AA^ effective interactions are renormalized 
by fitting the known binding energies of subsystems phenomenologically. Our analysis is 
performed systematically for ground and bound excited states of the series of a + x + A + A 
systems with no more adjustable parameters in this stage, so that these predictions offer 
an important guidance to interpret coming double-A events in the experiments and then to 
determine the level structure and the AA interaction unambiguously. 

In Section II, the calculational method with microscopic a + x + A + A four-body model is 
described. In Section III, the interactions are introduced. Calculated results are presented 
and discussed in Section IV. Summary is given in Section V. 
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II. MODEL AND METHOD 



In Ref. [|T3[ , the present authors already studied aaHg and I'jjie with the use of a + A + A 
three-body model and a + a + A + A four-body model, respectively. In the same manner, we 
study in this work the double-A hypernuclei /Jjie, ^1^1, aaLi, AA^i, aa-^^ AA^e on the 
basis of the a + x + A + A four-body model with x = n,p, d, t,^He and a, respectively. The 
d, t('^He) and a clusters are assumed to be inert having the (Os)^, (Os)^ and (Os)^ shell-model 
configurations, respectively, and are denoted by $s(a;) with spin s (= 1, | or 0, respectively). 

All the nine sets of the Jacobian coordinates of the four-body systems are illustrated in 
Fig. 1 in which we further take into account the antisymmetrization between two A particles 
and the symmetrization between two a clusters when x = a. The total Hamiltonian and 
the Schrodinger equation are given by 



H = T+Y.Vat + Kpauli , 
(a,fe) 



(2.1) 
(2.2) 



where T is the kinetic-energy operator and Vab is the interaction between the constituent 
particle-pair a and b. The Pauli principle between the nucleons belonging to a and x clusters 
is taken into account by the Pauli projection operator Vpauu which is explained in the next 
section as well as Vat- The total wave function is described as a sum of amplitudes of the 
rearrangement channels (c = 1 — 9) of Fig. 1 in the LS coupling scheme: 



^JM ( ^aZ) 



X AaSc 



EE E E 

c=ln,N,ul,L,XS,T,,I,K 



$,(x) [xi(Ai)Xi(A2)]i 



K 



(2.3) 



J JM 



Here the operator ^a stands for antisymmetrization between the two A particles, and 
Sa is the symmetrization operator for exchange between a clusters when x = a. xi{-^i) 
is the spin function of the i-th A particle. Following the Jacobian-coordinate coupled- 



rearrangement-channel Gaussian-basis variational method of Refs. [T^, |T8|, [T9|, ^ 



we take the functional form of 0„im(r), V'tvla/ (R) and ^^A/i(Pc) as 



= p"e-(''/'"')V,,(p), (2.4) 
where the Gaussian range parameters are chosen to he in geometrical progressions: 

Tn = ria""^ (n = 1 - rimax) , 

= R.A''-' (iV=l-iV^aJ, 

Pu = Pia"'^ (z/= 1 - Z/max) • (2.5) 

These basis functions have been verified to be suited for describing both short-range corre- 
lations and long-range tail behaviour of few-body systems |jl5|, [1^, 0, |18|, |19|, [2^, The 



eigenenergy E in Eq.(2.2) and the coefficients C in Eq.(2.3) are to be determined by the 
Rayleigh-Litz variational method. 

For the angular-momentum space of the wave function, the approximation with /, L, A < 2 
was found to be enough in getting satisfactory convergence of the binding energies of the 
states concerned presently. Note that no truncation is taken of the interactions in the 
angular-momentum space. As for the numbers of the Gaussian basis, n^^.^, A^max and I'rna.x, 
4 — 10 are enough. 

As far as the single A hypernuclei ^He, ^Li, ^Li, ^Li, ^Be and ^Be are concerned, the 
wave functions are described by Eq.(2.3) but with one of A particles omitted. As for the 
core nucleus itself, a + x, the wave function is given by 

^jM{a + x) =Y,CniSa^a)[^s{x)(j)ni{r)]jM ■ (2.6) 

n,l 

III. INTERACTIONS 

In the study of double-A hypernuclei based on the a -|- x -|- A -|- A four-body model, it 
is absolutely necessary and important to examine, before the four-body calculation, that 
the model with the interactions adopted is able to reproduce reasonably well the following 
observed quantities: (i) energies of the low-lying states and scattering phase shifts of the 
a + X nuclear systems, (ii) B\ of hypernuclei composed of x -|- A, x being d, t,^He, a, (iii) 
i?A of hypernuclei composed of a -|- x -|- A, x being n,p,d,t,^B.e, a and (iv) i?AA of \%Re 
= a-|-A+A. We emphasize that this severe examination were successfully done in the present 
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model as mentioned below. This encourages us to perform the four-body calculations, with 
no adjustable parameters at this stage, with high reliability of the results expected. 

A. Pauli principle between a and x clusters 

The Pauli principle between nucleons belonging to a and x clusters is taken into account 
by the orthogonality condition model (OCM)||22|. The OCM projection operator Vpauu is 
represented by 

yp,,H = lim A Y: l'^/(r«.))(</'/(rL)l (3.1) 

which rules out the amplitude of the Pauli-forbidden a — x relative states (j)f{rax) from 
the four-body total wave function |2^. The forbidden states are / = OS" for x = n(p), 
/ = {0^,0P} for X = rf, / = {0^, 1^, OP, OD} for x = t{^Re) and / = {0S,1S,0D} for 
X = a. The Gaussian range parameter b of the single-particle Os orbit in the a particle is 
taken to be 6 = 1.358 fm so as to reproduce the size of the a particle. The same size is 
assumed for clusters x = d,t and ^He to manage the Pauli principle avoiding the calculational 
difficulty. In the actual calculations, the strength A for Vpauu is taken to be 10^ MeV which 
is large enough to push up away the unphysical forbidden states in the very high energy 
region with keeping the physical states unchanged. Usefulness of this Pauli operator method 
of OCM has been verified in many cluster-model calculations. 

In some calculations |^, |IU], |rT|, |T^, ^ of three-body systems including two or three a 
clusters, use is made of an aa potential with a strong repulsive core so as to describe 
the Pauli exclusion role which prevents the two a cluster from overlapping. But, it is well 
known that this approximate prescription of the Pauli principle is not suited for the case 
where the presence of the third particle makes the two a clusters come closer to each other; 
in other words, the off-energy- shell behaviour of the repulsive potential is not appropriate in 
the three-body system. Moreover, there is no available potential reported for the a x systems 
{x = n,p,d,t and '^He) of this type. Therefore, we do not employ this prescription in the 
present systematic study of the structure change of the ax systems due to the addition of 
A particles. We take the orthogonality condition model instead which is suited even for the 
case of heavy overlapping between the two clusters. 



8 



B. ax interactions 



As for the potentials Vax between the clusters a and x, we employ those which have 
been often used in the OCM-based cluster-model study of light nuclei. Namely, they are 
the VaN potential introduced in Ref. the Vad and Vat potentials given in Ref.|T^ and 



the Vaa potential used in Ref. |g8| which reproduce reasonably well the low-lying states and 
low-energy scattering phase shifts of the a x systems. The potentials are described in the 
following parity-dependent form with the central and spin-orbit terms: 

i' 

•-max "max 

i=l i=l 
i" i'" 

max max 

i=l i=l 

(3.2) 

where / is the relative angular momentum between a and x, and s^^ is the spin of x. In 
the aa system the spin-orbit term is missing and the odd wave is forbidden by the Pauli 
principle. The additional Coulomb potentials are constructed by folding the pp Coulomb 
force into the proton densities of the a and x clusters. The parameters in Eq.(3.2) are 
listed in Table I (we slightly modified the strength of the central force in V^d and that of 
the spin-orbit force in V^t to obtain better agreement with the energy levels of ^Li and '^Li, 
respectively) . 

C. Ax interactions 

We derive the interaction between the A particle and the x cluster by folding the G- 
matrix type YN interaction (the YNG interaction) into the density of the x cluster in the 
same manner of our previous work on the double- A hypernuclei |T^ . The YNG interactions 
between A and are derived from the YN OBE models as follows: First the G-matrix 
equation is solved in nuclear matter at each kp, where the so called QTQ prescription is 
adopted for simplicity. Next the resulting G-matrix is simulated by a three-range Gaussian 



form with the strengths as a function of kp. Obtained YNG interactions are given in Ref. |29 

as 

1 + Pr 

[l^6,even + ^aa,even^ A ' N ) 

i=l 



6 

VAN {r; kp) = E [(4!Ln + ^aievenO-A " (T n) 



+ (^s:L + ^itdd'^A ■ ^n)^-^] e"""'' (3-3) 

where Pr is the space exchange (Majorana) operator. The strengths f are represented as 
quadratic functions of ; see Eq.(2.7) of Ref. and Table V of Ref. [|l^ for various 
original YN interactions. In the present work, we employ the Nijmegen model D interaction 
(ND). 

The Ax interaction is derived by folding the above v\N{f^'-, kp) interaction into the x- 
cluster wave function. The k-p depends on the mass number of the cluster x. Because of the 
operator Pj. in Eq.(3.3), the resultant Ax potential becomes nonlocal, the explicit form of 
which is given in Appendix of Ref. [|I^ . We summarize the functional form of the local and 
nonlocal parts of the A x potentials as 

yA.(r, r') = j^iV, + V/SA ■ s^) e-f^'^'dir - r') 

i=l 

+ E(^i + Uts^ ■ S,) e-7.(r+r')^-^.{r-r')^ ^ (3_4) 
i=l 

where sa = crA/2. Table II lists the parameters in Eq.(3.4) for (a) A a interaction, (b) 
A t(A^He) interaction and (c) A d interaction. They were determined in the following manner: 

i) Aa interaction: The AA^ spin-spin part vanishes by the folding into the a particle. The 
odd-force contribution is negligible to the A-binding energy of aHc. We determined the k^ 
parameter as k-p = 0.925 fm~^ in order to reproduce this binding energy (3.12 MeV) within 
the a + A two-body model. The AA^ odd- force having the same kp was determined by tuning 
the magnitude of foodd so as to reproduce, within the a + a + K model, the A-binding energy 
of the 1/2+ ground state of A^e. 

ii) Kd interaction: We determined the value of k-p = 0.84 fm~^ by fitting the experimental 
A-binding energy of the 1/2"'' ground state of within the d + A model where the AN 
odd force plays a negligible role. The odd force was determined, with the same kp kept, by 
reproducing the A-binding energies of the l/2f and 3/2f states of ^Li within the a + d + A 
model; we tuned v^^l^^ and wS^^dd- 

iii) At interaction: The experimental A-binding energies of the O"*" and l"*" states of 
were used to determine the even force of the AA^ interaction. The magnitude of kp and 
v^} even Were adjusted to reproduce the energies, kp being 0.84 fm~^. This value of kp was 
substituted into the kp used in the odd force of the AA^ interaction of the Ad interaction 
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with no other change. The resulting At interaction reproduces, by chance, the A-binding 
energy of the 1"^ ground state of \Li within the a + t + A model; the calculated energy is 
6.80 MeV while the observed one is 6.80 ± 0.03 MeV. 



D. AA^ interaction in jJ/^He (^^Li) 

In the study of aaHc (AA^i) with the a + A^ + A + A model and of the subsystem AHe(ALi) 
with the a + A^ + A model, it is inadequate to use the G-matrix type AA^ interaction because 
AA^ correlations are fully taken into account in our model space. Here, we employ a simple 
free-space AA^ interaction with a three-range Gaussian form, which simulates the Nijmegen 
model F (NF) AA^ interaction. Here, the SA^ channel coupling contribution is renormalized 
into the AA^ single channel with the closure approximation. The even- and odd-state parts 
of our AA^ interaction are represented as follows: 

1+Pr 



1=1 

1 p 

+ {Vi +Vi ' cta-ctn) — - — e . (3.5) 



First, the parameters are dertemined so as to simulate the AA^ scattering phase shifts calcu- 
lated with NF. Next, the second-range strengths and fg™'^''^ are adjusted so as to repro- 
duce the A binding energies of the 0^ and 1^ states of aH with the use of the A^ -|- A^ -|- A^ -|- A 
four-body model. Furthermore, strengths fg*^'^ and ^2'^*^''^ are adjusted within the framework 
ofa;-|-A-|-?T,-|-p four-body model so as to reproduce the observed binding energies of the 
ground-state spin doublet, 1/2"'" and 3/2"*" of aLI. Our resulting parameters in Eq.(3.5) are 
listed in Table HI. We further found that the energy of the ground state of AHe (aLI) mea- 
sured from the aHc — A^ threshold can be well reproduced with our AA^ interaction in the 
a + A^ + A three-body calculation; for ^He (aLI), the calculated energy is —0.17 MeV (0.57 
MeV), while the observed one is -0.17 MeV (0.59 MeV). 

E. A A interactions 

In the present model, since the AA relative motion is solved rigorously including the 
short-range correlations, it is not adequate to use the AA G-matrix interaction given in 



Ref. However, our AA interaction to be used in the present calculation should be 
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still considered as an effective interaction, since the couplings to SN and EE channels are 
not treated explicitly. Thus we employ the AA interaction represented in the following 
three-range Gaussian form: 

3 

VAAir) = J2{v, + <o-A ■ o-a) e-^''-'. (3.6) 

i=l 

It is enlightening here to keep some linkage to the OBE models in determining the inter- 
action parameters /i,, Vi and ff (z = 1 — 3). In our previous work on aaHc and \ jjie |jl3|, the 
interaction parameters were chosen so as to simulate the AA sector of the ND interaction 
which is a reasonable model for the strong attraction suggested by the old interpretation for 
double-A hypernuclei. The characteristic feature of ND is that there is only a scalar singlet 
instead of a scalar nonet, which gives strongly attractive contribution in AA as well as A^A^. 

The other versions of the Nijmegen models with a scalar nonet lead to much weaker 
AA attractions, which seems to be appropriate for the weak AA binding indicated by the 
NAGARA event. The NF is the simplest among these versions, which is adopted here as a 
guidance to construct our AA interaction: The outer two components of the above Gaussian 
potential {i = 1, 2) are determined so as to simulate the AA sector of NF, and then the 
strength of the core part {i = 3) is adjusted so as to reproduce the experimental value of 
i?AA(AAHe). The obtained values of parameters are given in Table IV. It is interesting that 
the resulting AA interaction is almost equal to the interaction obtained by multiplying a 
factor 0.5 on the above ND-simulating interaction employed in Ref. [O]. 



IV. RESULTS AND DISCUSSIONS 



Let us show the calculated results for a series of double-A hypernuclei with a + x + A + A 
structures (x = 0, n, p, ci, t,^ He, a) studied in the microscopic four-body cluster model. In 
order to understand the role of two A particles attached to the core nuclei, it is useful to 
compare the obtained level structures of the a + x + A + A double-A hypernuclei with those 
of the a + X nuclei and the a -|- x -|- A single-A hypernuclei. Then, we can see clearly how 
the ground and excited states of a -|- x nuclei are changed due to the participation of A 
particles. It should be noted again here that, in the model description of a + x + A + A, the 
observed low-energy properties of the a + x nuclei and the existing A-binding energies of the 
X + A and a + x + A hypernuclei have been reproduced accurately enough to give reliable 
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predictions for the double-A hypernuclei with no adjustable parameters of the interactions in 
the four-body calculations. It offers the most reliable ground for our cluster model that the 
threshold energies for every partition into sub-cluster systems are assured to be reproduced. 

A. Energy spectra 

In Figs. 2 — 7, the calculated level structure of a + x core nuclei, a + x + A hypernuclei 
and a + a; + A + A hypernuclei are illustrated side by side. There are shown all the ground 
and bound excited states of double-A hypernuclei predicted in the present model. In these 
figures, one sees clearly that injection of one and two A particles leads to stronger binding of 
the whole system and prediction of more bound states. But, there is no bound 'p-orbit' of A 
particle in single- and double-A hypernuclei with A < 10. In the bound states of double-A 
hypernuclei, two A particles are coupled to 5* = and therefore the spins and parities are 
the same as those of its nuclear core. 

Table V summarizes the calculated ground-state energies for the double-A hypernuclei 
including the 2+ excited state of aaBc. The results are expressed in terms of two quantities: 
One is the total energy measured from the breakup threshold ofa; + a; + A + A which is 
denoted as -Eaa- The other is i?AA which is the binding energy of two A particles with 
respect to the ground-state nuclear core a + x. 

The calculated values of -Baa can be compared with some experimental data, though 
they are quite limited at present. The most recent and clear data of the NAGARA event 
is used as a basic input of our model so that our AA interaction is adjusted to reproduce 
the experimental value B^^{j^jjie)= 7.25 ± 0.19±q;}^ MeV [||. It is of particular interest 
to compare the present result with another data which is not used in the fitting procedure. 
There is an event found in the E373 experiment, named as Demachi-Yanagi event |^0|, |3T| ; 
the most probable interpretation of this event is a production of a bound state of \/Jie 
having I^aa' = 12.33±Q;|f MeV which is obtained by assuming = 0.15±q^ MeV. In the 
emulsion analysis there is no direct evidence for the production of aa-^^ excited state. 
However, if the produced AA^e is interpreted to be in the ground state, the resultant AA bond 
energy becomes repulsive, contradictorily to the NAGARA event. From the viewpoint of the 
present study, the Demachi- Yanagi event can be interpreted most probably as observation 
of the 2+ excited state in aa-^^j calculated value of B {^!yQ^jJ^e{2^) ) is 12.28 MeV, which 
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agrees with the above experimental value. This good agreement suggests that our level 
structures calculated systematically are predictive and useful for coming events expected to 
be found in the further analysis for the E373 data. Now it should be stressed that the above 
experimental data of Be(2+) leads to no information for the ground-state value of i?AA 
unless the theoretical value (2.86 MeV in our case) of the excitation energy of ^°Be(2+) is 
utilized. 

On the other hand, the old experiment by Danysz et al. on the pionic decay of 
jv°ABe(0+) ^ \Be{l/2+) + j9 + vr" gave 5f/(\o^Be(0+)) = 17.7 ± 0.4 MeV. This value has 
been used for a long time, which means the strongly attractive AA interaction. However, it 
should be noted that the authors also suggested the possibility of another decay ^^36(0+) 
iBe(3/2+,5/2+) + p + vr^ (Table 5 of Ref.0); the same was pointed out in Ref. 0, too. 
In this case, the value of -B^a (aaBb (0^)) is modified to 14.6 ± 0.4 MeV, which is obtained 
by using the excitation energy of ABe(3/2+, 5/2+) to be 3.05 MeV ||32|. This modified 
value turns out to be not contradictory with our calculated value, 15.14 MeV. A similar 
re-interpretation, with the hypernuclear excited states taken into account, may be needed 
also for the E176 event which was identified as ^a^ (aa^g) with the strongly attractive 
(repulsive) AA interaction. 

Thus, we have understood the consistency between the experimental data and our theo- 
retical results of aa-^s. We, therefore, discuss on the level structures of double-A hypernuclei 
in more detail. As seen in Figs. 2 — 7 and Table V, the A particle plays a glue-like role so 
that a whole system becomes of stronger binding. This effect in a double-A nucleus is more 
enhanced than that in the corresponding single-A nucleus. One can see a typical example in 
the case of AA-^i in Fig- 3. For the unbound nuclear system of ^Li, a single A cannot make 
a bound system of aLI, but, the addition of one more A particle leads to a bound system of 
^^Li whose ground state is of weakly-binding with respect to the ^pYie+p threshold. 

The bound excited states of double-A hypernuclei predicted in the present cluster model 
are summarized as follows: In aaHs and AA-Li, the ground states are both bound but no 
excited states are predicted. It is needless to say, there are no bound excited states in 
double-A hypernuclei with A <Q since there is no bound excited state in their core nuclei. 
The lightest double-A hypernucleus that has at least one excited state is AA^i- In aa^^ 
we predict two T = excited states in the bound-state region. It is expected to have a 
T = 1,0+ bound excited state in AA^i which corresponds to the T = 1,0+ state in ^Li 
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nucleus at = 3.56 MeV, but the state is not shown in Fig. 4 because the T = 1 state 
may have five-body structure and is out of scope of the present cluster model. We predict 
three bound excited states in aaLI (yviBe). There is one bound excited state in ^A^e as 
mentioned before. It will be challenging to discover these excited states one by one as well 
as the ground states. 



B. Dynamical change of the core nucleus 

It is interesting to look at the dynamical change of the a + x nuclear cores which occurs 
due to the successive participation of two A particles. The possibility that a nuclear core 
shrinks due to an attached A particle has been theoretically pointed out using the a + x + A 
cluster model of light p-shell A hypernuclei [^. An updated prediction |Q was made 



specifically on a shrinkage in size by 21% to be seen in ^Li. The recent measurement of 
7-ray transition rate in ^Li has confirmed quantitatively the shrinkage effect predicted 



in both the old calculation and the updated one. It is quite reasonable, therefore, that 
in a double-A hypernucleus the participation of one more A particle can induce further 
shrinkage of the nuclear core. Such an effect has been also investigated systematically using 
the molecular orbital model for ^A^Be {n = l — 4:) = a + a + nA ^6 . 



In order to see such shrinkage effect, we show three physical quantities: first in Table 
VI we list the r.m.s. distance between a and x, fax- As the number of the A particles 
increases, fax turns out to shrink significantly due to the glue-like role of the attached A 
particles. For example, one sees fax changes as 4.11 —>■ 3.44 —>■ 3.16 fm for ®Li —>-\Li^ 
aaLI. Participation of the second A gives rise to about 8 % reduction of fax except x = n. 
Secondly, in more detail, it is worthwhile to demonstrate in Fig. 8 the change of the a — n 
two-body density (correlation function) p{ran) in the ^He, and AA^e when A particles 
participate successively, which again manifests the shrinkage effect. Thirdly, this shrinkage 
effect is seen in the large change of the expectation value of the relative kinetic energy, 

< Tax >, and that of the potential energy, < Vax >, in the a — x subsystems. When the 
a and x clusters approach to each other, the increase of < Tax > overcomes the gain of 

< Vax >, and the sum < Tax + Vax > increases appreciably. In spite of this energy loss in 
the a — x core system, the core shrinkage is realized by the stronger energy gain of the A — a 
and A — x parts. 



15 



C. AA bond energy 



In Fig. 9 we reveal the contributions of the AA interaction to the total binding energies 
of double- A hypernuclei ^^^Z. Here the calculated values of Saa(aaZ) in the ground states 
are shown by closed circles. In order to extract the contribution of the A A interaction, we 
perform the same calculations by putting Vaa = 0. The obtained values are denoted as 
-^aa(aaZ; ^aa = 0) and shown by open circles in the figure. It should be noted that the 
effect of the dynamical change of the a -\- x core due to the AA^ interactions is included 
in the four-body estimate of Saa and -Baa (Vaa = 0). Since the AA interaction is not so 
strong compared with the AA^ interaction, the core-rearrangement effects included in Saa 
and -Baa(14.a = 0) are similar to each other. Then, naturally the pure effect of the AA 
interaction is given by the difference 

(^aZ) = Saa(^aZ) - Saa(^aZ; l^AA = 0). (4.1) 

We consider Vaa''^ as the AA bond energy which should be determined essentially by the 
strength of the AA interaction. Now in Fig. 9, we find that the magnitude of Vaa"^) the 
energy difference between the closed and open circles, is almost constant at ~ 1 MeV for all 
the double- A hypernuclei with A = 6 — 10. The detailed values of Vaa^^ are fisted in Table 
V. 

So far the following intuitive formula has been often used to estimate the AA interaction 
strength: 

ABaa(^aZ) ^ i?AA(^AZ) - 2Sa(^-'Z). (4.2) 

It is worthwhile to point out the problems underlying in this formula: This expression 
includes three problems which come from i) the mass-polarization term of the three-body 
kinetic-energy operator, ii) the AA^ spin-spin interaction and iii) the dynamical change of 
the core nuclear structure. 

The problem i) is stated as follows: In the a -|- A -|- A three-body model for aaHc (gener- 
ally, "a" may be replaced by "spinless frozen-core nucleus"), if one takes the non- Jacobian 
coordinate set FqAi and r^Aj, the Shrodinger equation may be written, in a self-explanatory 
notation, as 

[ --^ V„Ai - TT- Vo,A2 - — VaAi ' VaAi 
-^PaAi '^paA2 '"'a 

+ K.A1+KA2 + Vk,k, - E \^jM{l^^e) = 0. (4.3) 
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If the third term of the kinetic energy, the so-called mass-polarization term, and VaiA2 
are neglected, we have the trivial solution ~E{= -Baa) = 2B\. Therefore, the quantity 
Ai?AA = i?AA — 2i?A stands for the contribution from the neglected two terms. In aa^^, 
the contribution to i?AA from the mass-polarization term is +0.13 MeV which explains the 
difference between ABaa = 1.01 MeV and the AA bond energy V^^"'^ = 0.88 MeV in Table 
V. This contribution decreases rapidly as the core-nuclear mass increases (+0.01 MeV in 
aW- 

Next, we discuss about the second problem, an effect of the AA^ spin-spin interaction on 
ABaa of Eq.(4.2). In Fig. 10, the calculated values of ABaa are illustrated by the dashed 
bars. One notices clearly that Ai?AA has peculiar mass dependence in which some interesting 
mechanism is included. It should be remarked here, however, that, as was already pointed 



out by Danysz et al. 0, the traditional definition of Eq.(4^) is of simple meaning only when 
the nuclear core is spinless. On the other hand, in the case of nuclear core with spin, the 
single-A binding energy B\ to be subtracted from i?AA is distributed over the ground-state 
doublet of the corresponding single-A hypernucleus. 

Here, we remark the fact that the AA^ spin-spin interaction is not effective (cancelled 
out) in the double-A hypernuclei having the AA spin-singlet pairs. In the parent single-A 
hypernuclei, however, the spin-spin interaction plays an important role in giving rise to 
the energy splitting of the ground-state doublet. The typical and unique example known 
experimentally is the spin-doublet in ^Li with J = 1^ (ground; B^ = 5.58 MeV) and J = 
^~^{Ex = 0.69 MeV; i?A = 4.49 MeV). Considering this effect, one should use the spin- 
averaged value i?A(ALi) = |-SA(|gg ) + |-Ba(|^) instead of -BA(|gg ) when one likes to deduce 
A_Baa from the AALi(l''') ground state data, if any. If we adopt this prescription also for the 
adjacent systems, we may use 



Ba 


:iHe) = 


i^A(lg.s. 


) + IBa{2-) 


= 4.02 MeV, 


Ba 


(iLi) = 


\Ba{11.) 


+ \Ba{2-) 


= 4.31 MeV, 


Ba 


(iLi) = 


I^a(|^,.) 


+ iBAit) 


= 5.12 MeV, 


Ba 


(iLi) = 


i^A(l^s.) 


+ IBa{2-) 


= 6.58 MeV, 


Ba 


(iBe) = 


i^A(lg;s. 


) + !i?A(2-) 


= 6.48 MeV. 
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Here, Sa of the excited states are taken from our calculation. In general, we have 



+ 



2Jo + l 
Jo + 1 
2Jo + l 



-Ba(a ^Z; Ji = Jo - -) 
-Sa(a~^Z; Ji = Jo + -), 



where Ji = Jq^ \ denote the doublet spins of the q;+x+ A system, Jo being the ground-state 



spin of the a -\- x nuclear core. For the two spinless cases {x — Q and ct), needless to say. 



In Fig. 11, the sohd bars illustrate ASaa- Though ASaa is free from the effect of the 
AA^ spin-spin interaction, its magnitude for ^4 = 7 — 10 deviates significantly from ASaa 
(AAHe)= 1.01 MeV. The deviation comes from the effect of the dynamical change in the core 
nucleus structure (shrinkage in the a — x distance) due to the partition of the A hyperons, 
and turns out to be maximum in the case of aaBc. We emphasize that, even if one employs 
ASaa, it is impossible to extract any consistent value of the AA bond energy from Fig. 11 
in which ASaa scatters in a range of a factor of two. 

As mentioned above, a consistent estimation of the AA bond energy (0.9 — 1.0 MeV, 
nearly independent of the mass number, as seen in Table V) can be obtained by taking 
Vaa^'^ of Eq. (4.1) as the definition of that energy, though help of the theoretical calculation 
with Vaa = is necessary. 

V. SUMMARY 

We have carried out structure calculations of AA^e, AA^e, AA^i, AA^i, AA^i , AA^e and 
^^Be taking the framework of a-\-x-\-h.-\-h. model with a; = 0, n, p, d, t,^He and a, respectively. 
We determined the interactions between constituent particles so as to reproduce reasonably 
the observed low-energy properties of the a + x nuclei and the existing data of A-binding 
energies of the a; -|- A and a -\- x -\- A. systems. The A A interaction was constructed so as 
to reproduce the SAA(AAHe) given by the NAGARA event within our a -|- A -|- A model, 
where the long-range part of our interaction was adjusted to simulate the behavior of the 
appropriate OBE model (NF). With no adjustable parameters, the four-body calculations of 




ASaa(aaZ) = -Baa(aaZ) - 25a(a ^Z). 



(4.4) 
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the a+x+A+A systems were performed accurately using the Jacobian-coordinate Gaussian- 
basis coupled-rearrangement-channel method. Obtained energy spectra of the double-A 
hypernuclei with A = 6 — 10 are summarized in Fig. 12. 
Major results to be emphasized here are as follows: 

(1) It is striking that the calculated -Baa of the 2+ excited state in AA^e, 12.28 MeV, 
agrees with the experimental value -Baa (aaBg) = 12.33±o2f MeV in the Demachi-Yanagi 
event |^0|, We therefore interpret this event as observation of the 2"*" excited state of 



10 
AA 



Be. The agreement suggests that our systematic calculations are predictive for coming 



events expected to be found in the further analysis of the E373 data, etc. 

(2) Together with the energy spectrum of each double-A hypernucleus, those of the cor- 
responding core nucleus and single-A hypernucleus are exhibited side by side in Figs. 2 — 7 
so as to see clearly that injection of one and two A particles leads to stronger binding of 
the whole system and prediction of more bound states. In the bound states of any double-A 
hypernucleus, two A particles are dominantly coupled to S" = and hence the spin and 
parity become the same as those of its nuclear core, but the theoretical -Baa values are of 
importance to guide the analysis of the emulsion experiments. 

(3) Dynamical change of the a + x nuclear core by the participation of the A particles 
is substantially seen in double-A hypernuclei; there occurs, averagely speaking, about 8 % 
shrinkage of the a — x distance compared with the distance in the single-A hypernucleus. 
This shrinkage is realized by the large energy gain in the A — a and A — x parts which 
overcomes the energy loss in the a — x relative motion. 

(4) We estimated the AA bond energy using the faithful definition Vaa**^ = -Baa — 
5aa(Vaa = 0) and found it to be 0.88 MeV for ^^He and 0.93 - 0.98 MeV for the other 
double-A hypernuclei. We demonstrated that the quantity A-Baa = -Baa — 2-Ba is not a good 
measure of the AA bond energy since A-Baa is free from neither the contribution from the 
splitting of the ground-state doublet in the single-A hypernucleus nor that of the structure 
change of the core nucleus. In fact, the value of A-Baa scatters from 0.28 to 1.68 MeV for 
the double-A hypernuclei with A = 6 — 10. We then modified A-Baa by A-Baa = -Baa — 2-Ba 
with -Ba being the spin-average of -Ba's for the ground-state spin-doublet. We found, how- 
ever, that A-Baa still ranges from 0.83 to 1.68 MeV due to the structure change of the core 
nucleus. Direct use of -Baa itself rather than the use of A-Baa or A-Baa is recommended 
when the experimental result and calculational result are compared to each other. 
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In conclusion, the present precise and extensive four-body cluster-model calculation can 
be an opening of the spectroscopic study of double- A hypernuclei. 
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TABLE I. Parameters of (a) aa interaction, (b) at (A^He) interaction, (c) ad interaction 
and (d) aN interaction defined in Eq.(3.2). Size parameters are in fm"^ and strengths are 
in MeV. The ^-S'o scattering length is -0.575 fm and the effective range is 6.45 fm. 





(a) a a 


interaction 




I 


1 


2 


3 


A 


0.1111 


0.2777 


0.3309 


Vi 


-1.742 


-395.9 


299.4 




0.0 


0.0 


0.0 




(b) a t (ct^He) interaction 




i 


1 


2 


3 


A 


0.0913 


0.1644 


0.2009 




6.9 


-43.35 


-51.7 


Pf 


0.0913 


0.1644 


0.2009 


vr 


6.9 


43.35 


-51.7 


li 


0.28 






vt 


-1.2 








0.28 








1 9 
1 .z 








(c) a d 


interaction 




i 

A 
V 


1 

0.2 

-64.21 






P\ 


0.2 






vf 


-10.21 






li 


0.3 
-4.0 






n 

7f 


0.3 






y-ls,p 

^ i 


-4.0 








(d) aN 


interaction 




i 

A 

Vi 


1 

0.36 
-96.3 


2 

0.9 
77.0 


3 


Pf 


0.2 


0.53 


2.5 


vr 


34.0 


-85.0 


51.0 


li 


0.396 


0.52 


2.2 


vt 


-20.0 


-16.8 


20.0 


7f 


0.396 


2.2 






6.0 


-6.0 
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Table II. Parameters of (a) A a interaction, (b) Ai(A^He) interaction and (c) A d interaction 



defined in Eq.(3.4). Size parameters are in fm ^ and strengths are in MeV. 





(a) A a 


interaction 




i 


1 


2 


3 




0.2752 


0.4559 


6123 




-17.49 


-127.0 


497.8 




0.0 


0.0 


0.0 


It 


0.1808 


0.1808 


1 808 


Si 


0.4013 


0.9633 


2.930 


Ui 


-0.3706 


-12.94 


-331.2 




0.0 


0.0 


0.0 




(b) A t (A^He) interaction 




i 


1 


2 


3 


Hi 


0.2874 


0.4903 


6759 


V 


-14.16 


-108.0 


425.9 




2.379 


10.91 


-126.9 


li 


0.2033 


0.2033 


U.zUoo 




0.3383 


0.8234 






-0.2701 


-9.553 


-231.6 




-0.2615 


1.433 


97.05 




(c) Ad 


interaction 




i 


1 


2 


3 


A 


0.3153 


0.5773 


0.8532 




-10.84 


-88.36 


167.2 




2.734 


14.35 


-179.9 


li 


0.2710 


0.2710 


0.2710 


Si 


0.2470 


0.4870 


1.924 




-0.1862 


-5.844 


-3.065 


U! 


-0.2705 


1.566 


100.4 



TABLE III. Parameters of the A interaction defined in Eq.(3.5) which is used only in the 
q; + A^ + A and o; + A" + A + A systems {x = N). Size parameters are in fm~^ and strengths 
are in MeV. 





A N interaction when x — N 




i 


1 


2 


3 


l^i 


0.5487 


1.384 


6.250 


„ .even 


-10.40 


-87.05 


1031 


even,<T 
Vi 


0.2574 


17.09 


-256.3 


„,odd 


-5.816 


-18.29 


4029 


odd, IT 
Vi 


-0.959 


-9.184 


-573.8 
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TABLE IV. Parameters of the A A interaction defined in Eq.(3.6). Size parameters are in 



fm ^ and strengtlis are in MeV. The ^5*0 scattering length is -0.575 fm and the effective 
range is 6.45 fm. 



A A interaction 

i I 2 3 

/ii 0.555 1.656 8.163 

Vi -10.67 -93.51 4884 

< 0.0966 16.08 915.8 



Table V. Calculated energies of the ground states of A = 6 — 10 double-A hypernuclei based 
on the a + X + A + A four-body model (x = 0, ra, p, rf, t^^ He, and a) . E\\ are measured from 
the a + X + A + A threshold. The AA bond energy V^^'^'^ is defined by Eq.(4.1). Information 
on the 2+ excited state of jjj^e is specially added so as to demonstrate the agreement with 
the experimental result. 







Eaa 
(MeV) 


Baa 
(MeV) 


-^AA 

(MeV) 


1 ;bond 
I^AA 

(MeV) 


A^He 


0+ 


-7.25 


7.25 


7.25 ±0.19" 


0.88 


aIHc 


2 


-8.47 


9.36 




0.96 


aaLi 


3- 
2 


-7.48 


9.45 




0.95 


AA^l 


1 + 


-12.10 


11.44 




0.98 


AA^l 


3- 
2 


-17.05 


14.55 




0.98 


AA^6 


3- 
2 


-16.00 


14.40 




0.97 


10 -n 
AA^6 


0+ 


-15.05 


15.14 


17.7 ±0.4^ 
14.6 ±0.4^ 


0.93 


10 Rp 

AA^6 


2+ 


-12.19 


12.28 


12.33±[]:^f ^ 


0.93 



" Ref.0. 

Ref.PI. Also see text for the second value. 



Ref. 30, 31 



Table VI. Calculated r.m.s. distances between a and x, Tax, in core nuclei, single A hyper- 
nuclei and double-A hypernuclei (x = n,d,t,a). The expectation values of kinetic energy 
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and potential energy between a and x, < Tax >, < Vax > and < T^x + Vox > are also listed. 
For ^He and ^Be, fa-x are not calculated since they are resonant states. 











^ Tax ^~ ^ax ^ 


iHe 


5.79 


7.86 
11.38 


-6.97 
-9.92 


0.89 
1.46 


aIHc 


3.92 


15.19 


-11.95 


2.24 


^Li 


4.10 
3.44 


11.59 
15.59 


-13.06 
-16.70 


-1.47 
-1.11 




3.16 


18.86 


-19.54 


-0.68 


^Li 


3.69 
3.30 


17.45 
21.85 


-19.95 
-24.00 


-2.50 
-2.15 


AA^l 


3.05 


26.74 


-28.33 


-1.59 




3.78 


7.21 
14.90 


-7.12 
-14.14 


0.09 
0.76 


i^Be 


3.44 


19.49 


-17.96 


1.53 



Figure captions 

FIG. 1. Jacobian coordinates for all the rearrangement channels {c — 1 — 9) of the 
q; + x + A + A four-body system. Two A particles are to be antisymmetrized, and a and x 
are to be symmetrized when x — a. 

FIG. 2. Calculated energy levels of ^He, ^He and ^iHe on the basis of the a+n, a+n+A 
and a + n + A + A models, respectively. The level energies are measured from the particle 
break-up thresholds. 

FIG. 3. Calculated energy levels of ^Li, ^Li and aa^I on the basis of the a-f-p, a-|-p-|-A 
and a + p + K + K models, respectively. The level energies are measured from the particle 
break-up thresholds. 
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FIG. 4. Calculated energy levels of ^Li, ^Li and ^\Li on the basis of the a + d, a + d + A 
and a + d + A + A models, respectively. The level energies are measured from the particle 
break-up thresholds or are given by excitation energies Ey^. 

FIG. 5. Calculated energy levels of ^Li, |Li and ^^Li on the basis of the a + t, a + t + A 
and a + t + A + A models, respectively. The level energies are measured from the particle 
break-up thresholds or are given by excitation energies E^. 

FIG. 6. Calculated energy levels of ^Be, ^Be and ^A^e on the basis of the a + ^He, 
a + ^He -|- A and a + '^Hc + A -|- A models, respectively. The level energies are measured 
from the particle break-up thresholds or are given by excitation energies £^x- 

FIG. 7. Calculated energy levels of ®Be, ^Be and l\Be on the basis of the a+a, a+a+A 
and q; -I- q; -I- A -I- A models, respectively. The level energies are measured from the particle 
break-up thresholds or are given by excitation energies Ey^. 

FIG. 8. The a — n two-body densities (correlation functions), p{ran), of ^He(3/2~), 
^He(l~) and ^^He(3/2~). Here, it is multiplied by r^„. 

FIG. 9. Calculated values of Saa(^^Z) in the ground states are given by closed circles. 
The same quantities but calculated by putting Vaa = 0, namely Saa(aa2; Vaa = 0), are 
shown by open circles. 

FIG. 10. Calculated values of A5aa(aa^) defined in Eq. (4.2). 

FIG. 11. Calculated values of A5aa(aa^) defined in Eq. (4.3). 

FIG. 12. Summary of the energy levels of the double-A hypernuclei aaHc, AA^e, AA^i, 
AALi, AA^i, AA^e and AA^e calculated using the a + x + A+ A model with a; = 0, n,p, d, t,^He 
and a, respectively. 
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